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Abstract 

Future lensing surveys will be nearly full-sky and reach an unprecedented depth, probing scales closer and 
closer to the Hubble radius. This motivates the study of the cosmic shear beyond the small-angle approxi- 
mation and including general relativistic corrections that are usually suppressed on sub-Hubble scales. The 
complete expression of the reduced cosmic shear at second order including all relativistic effects was derived 
in [T]. In the present paper we compute the resulting cosmic shear bispectrum when all these effects are 
properly taken into account and we compare it to primordial non-Gaussianity of the local type. The new gen- 
eral relativistic effects are generically smaller than the standard non-linear couplings. However, their relative 
importance increases at small multipoles and for small redshifts of the sources. The dominant effect among 
these non standard corrections is due to the inhomogeneity of the source redshift. In the squeezed limit, its 
amplitude can become of the order of the standard couplings when the redshift of the sources is below 0.5. 
Moreover, while the standard non-linear couplings depend on the angle between the short and long mode, the 
relativistic corrections do not and overlap almost totally with local type non-Gaussianity. We find that they 
can contaminate the search for a primordial local signal by ^ 10- 



1 Introduction 

With the advent of future wide-field surveys, weak gravitational lensing will become a premier probe of 
cosmology and an important tool to constrain dark energy, neutrinos and the initial conditions (see for 
instance |U In order to fully exploit the potentiality of the convergence and shear fields, it will be 
important to use their whole statistics. In particular, the lensing bispectrum represents a complementary 
probe to the power spectrum, as it will provide constraints that are comparable to those obtained from the 
power spectrum alone [1]. 



One of the primary interests of wide-field surveys is to look for primordial non-Gaussianities generated in 
the early universe and the lensing bispectrum represents a natural tool to capture such features [3 |6l [71 HI |9] . 
Given the current and forecasted constraints on /nl, there is a plethora of second-order effects intervening 
between the initial conditions and the observations that could be potentially relevant. Some of these effects, 
such as the second-order growth of matter fluctuations in Eulerian perturbation theory, the Born correction, 
the lens-lens coupling and the non-linear conversion between the galaxy shape distortion and the (observable) 
reduced shear, have been thoroughly studied in the past years and are expected to dominate the lensing 
3-point statistics on small-angular scales piIl[TT|P^[T51[n i. [T5 1. [ro i [T7 1 [T5 1 [Tg]. 

When studying weak lensing one usually restricts the analysis to small angular separations (or large 
multipole moments), which is justified by the fact that so far cosmic shear surveys have covered only a limited 
portion of the sky. However, nearly full-sky surveys are currently under preparation and on large angular 
scales general relativistic second-order couplings will become relevant. These corrections are usually neglected 
because they are suppressed by the ratio between the scale probed and the Hubble scale, but are important 
on angular scales of the order of the angular diameter distance to the source. These are analogous to the 
second-order effects from general relativity affecting the CMB bispectrum when at least one of the scales 
probed is larger than the horizon at recombination PUIIHIIH]. In the CMB, these effects have been recently 
found to be negligible for the contamination of Planck searches for a primordial local signal in the squeezed 
limit [23 In case of cosmic shear observations, these effects however cannot be a priori neglected. The 
aim here is to provide an exhaustive evaluation of the impact of these contributions to observations. Those 
results extend recent studies on the impact of relativistic corrections to the observations of large-scale galaxy 
clustering in [351 [Ml [571 [HI HHj . 

The complete study of all second-order effects in the cosmic shear, including the general relativistic ones, 
has been recently undertaken in IJ by solving at second order the Sachs equation [30 , which describes the 
deformation of the cross-section of a light bundle, mapping galaxy shapes into their angular images. The 
main results of this paper are reviewed in Sect.[2l 

The goal of the current paper is to compute the bispectrum of the electric part of the cosmic shear from 
all second-order effects and compare the importance of the general relativistic corrections with the standard 
couplings. We review the computation of the shear power spectrum in Sect. [3] where we also present one of the 
numerical shortcuts we will use throughout the paper, the second-order Limber approximation. The formal 
expression of the bispectrum, computed in Sect. [4j is determined by two types of non-linear contributions: 
the dynamical couplings, which depend on the particular metric solution of Einstein's equations at second 
order, and the geometrical couplings, so-called because they depend on the geometry of the solution of 
the Sachs equation at second order. Among the latter we also include the corrections coming from the 
inhomogeneities of a fixed-redshift source plane. The calculations of Sect. [4] involve multiple integrals and 
complicated manipulations of spherical harmonics and Wigner symbols. (In appendix [C] one can find useful 
relations for these calculations.) We are afraid to say that those calculations are extremely lengthy and 
technical. Impatient readers can skip this section and go directly to Sect. ^ In this section we present our 
results concentrating on the squeezed limit, where one of the multipole moments I involved in the bispectrum 
is small and we discuss the functional forms of the resulting terms as well as their relative importance. In 
Sect. |6] we compare our results with the bispectrum generated by primordial non-Gaussianities of the local 
type and we compute their contamination to a primordial /^l. Finally, in Sect, [t] we conclude and discuss 
the results of the paper. 
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2 Full-sky lensing shear at second order 



2.1 The shear and spin-s spherical harmonics 

The shear is characterized by a (2 x 2)-traceless and symmetric tensor 7^^, whose components are defined with 
respect to a particular choice of axes about each direction on the sky. This tensor describes the deformation of 
the image of distant galaxies by the gravitational field of the cosmic structures intervening between emission 
and observation. The two independent components of jab, 7i and 72, can be conveniently combined into a 
single complex field representing the shear, 

7 = 71 + «72 • (1) 

Under a right-handed rotation of the axes by an angle a about the line of sight h, this complex field transforms 
as 7 — > e*^"7. Thus, it represents a spin-2 field that can be appropriately expanded in terms of spin-weighted 
spherical harmonics, as 

7("-) = X! 2"'™ 2Yijn{n) ■ (2) 

Im 

The complex conjugate of 7, 7* = 71 — 172, is then a spin-(— 2) field, and can therefore be expanded as 

7*("-) = X! -20;™ -2Yim{n) ■ (3) 

Im 

One can introduce spin raising and lowering operators that can be used to relate quantities of difi^erent 
spin. The spin raising operator is denoted as and the spin lowering operator as 0. For a general spin-s 



field sX, these are defined as 



sX = - sin' 9{dg + i^^^d^ ] (sin-^ 6) ,X , (4) 



smf 



0,X = - sin-' eidg- i^d^ ] {sin' 9) sX . (5) 
' sm" ' 



Note that and commute only when applied to a spin-0 quantity. In general one has 

{00- 00),X ^2s,X . (6) 

These operators can be used to obtain spin-0 quantities. Acting twice with and respectively on 7 and 7* 
in eqs. ^ and ([s]) , and using the orthogonality properties of the spherical harmonics one obtains an expression 
for the expansion coefficients 2aim and _2a/m as 

2ai„, = l{l, -2) j dh Yi^Jn) 0^j{n) , (7) 

-2aim - l{l,~2) J dn Yi*^{n) 0^j*{h) , (8) 
where for convenience we have defined I as 



As for the Cosmic Microwave Background (CMB) polarization, we introduce parity eigenstates, spin-0 
quantities, called the "electric" and "magnetic" parts of the shear, defined as [5T1 [52] 

E{n)=J2aE,hnY,m{h) , (10) 
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B{h) = '^aB,lrn Yirn{h) , (11) 
Im 

where the coefRcients aE,im and as^im are given by 

aEJm = —-^{-2aim + 2a/m) , (12) 
1 

a-Bdm = —-^{-2aim — 2a/m) ■ (13) 

Under parity transformation, E and B change as E E and B —B. Furthermore, in contrast to the shear 
components 71 and 72, these two spin-0 quantities have the advantage of being rotationally invariant. 

2.2 The shear up to second order 

We consider a flat FLRW metric background. In [1] we have computed the shear at second order by using 
a perturbed metric in the so-called generalized Poisson gauge. In this gauge the second-order metric reads 

[33 [S] 

ds^ = a^irj) [-e^'^dri^ + 2w, d-qdx' + (e^^'^J.j -f h,j) dx' dx^] , (14) 

where the vector component uoi is divergenceless, diUJi — 0, and the tensor component hij is divergenceless 
and traceless, dihij = = ha. While the scalar perturbations 4> and "0 contain first- and second-order 
contributions, since we neglect primordial vector and tensor perturbations, uii and hij are only second-order 
quantities. 

Note that we have chosen to write the gravitational potentials in the metric in the exponential form. This 



choice is convenient for two reasons. First, up to second order, in this form the metric (14) is conformal to 
ds^ = —e^^'^'^^^drf + 2uii drj dx^ + {5ij + hij) dx^ dx^ so that the effect of scalar perturbations on null geodesies 
will be only through the Weyl potential, defined as 

* = (0 + V)/2. (15) 

The second reason, which will be explained in more details in section |2.2.2[ is that with this choice the 
relativistic second-order contributions to cf) and ?/; vanish in the squeezed limit. 

In this metric the Icnsing shear is obtained by solving Sachs equation |30| , which describes linear deforma- 
tions of the infinitesimal cross-section of a light bundle in the optical limit, and maps galaxy intrinsic shapes 
into their angular images. At first order one finds the full-sky expression of the shear fielc|^ 

l{n)= r dx^^^$^^{x.x) , (16) 
Jo XsX 

where we have conveniently defined % = 770 — where 770 is the conformal time today, so that x = fix is the 
background photon geodesic. The subscript S in Xs refers to the source. 
At second order, there are four sources of non-linearities [Tj: 

1) The mapping solution of the Sachs equation is linear in the angular deformation, but it is non-linear in 
the Weyl potential encountered by the photon from emission to observation. This induces non-linearities 
even when using the linear part of the metric. 



^Surprisingly, to our knowledge this full-sky expression was explicitly derived only recently in 1 although it is 
implicit in [35] and in the context of CMB lensing, see for instance [361 137) . The form was painstakingly rederived in 
[38) in the context of full-sky cosmic shear observations. 
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2) What we observe is actually the reduced shear [THl [IHl HZ] j i-e. the ratio between the anisotropic and 
isotropic deformations, and there are non-linear corrections introduced when taking this ratio. These 
will be proportional to the product of the first-order shear and the first-order convergence. 

3) Observationally, we are mapping galaxies located at a given redshift zs- Thus, we expect second-order 



contributions to the observed shear coming from perturbing zs in the first-order expression ( 16 ) 



4) The metric contains second-order terms in the initial conditions. This induces non-linearities even when 
using the linear solution of the Sachs equation. 

The first three contributions are independent of the second-order components of the metric: we will 
collect them under the name of "geometrical" contribution, because they depend on the geometry of the 
linear mapping, solution of the Sachs equation at second order. On the other hand, the fourth contribution 
depends only on the linear mapping at first order, but depends on the particular metric solution of the 
Einstein equations at second order. Thus, it is natural to call it the "dynamical" contribution. Note that this 
separation depends on the particular gauge chosen, hence it is not completely physical. However, for reasons 
that will appear more clearly in the following, it is a convenient distinction that we adopt hereafter. 

2.2.1 Geometrical couplings 

Except for the contribution coming from the perturbation of the redshift, which we discuss below, the geo- 
metrical couplings depend exclusively on the Weyl potential 5*. As explained above, this is a consequence 



of the choice of the exponential of and ^ in the form of the metric in eq. (14). Some of the geometrical 
couplings dominate in the large-^ limit - corresponding to small angular scales - because they contain more 
operators $ or $. Indeed, in harmonic space to each of such operators is associated an Z-factor in the angular 
power spectrum and bispectrum. These dominant terms, that are standard in the literature, can be written 
as 

I'gZ^ = - Prf^ r d^^^^^^^\$ ($H{x.x)$^{^ ,x') + $$^{x,S)$'^W,x')) 
Jo Jo XsX XX L V J 

+ 2^^'(x,af)^*(x',x')] (17) 

+ T'dx r dj^^^^^^^^$$^{x.x)$H{x:,x') . 
Jo Jo XsX XsX 

The first two lines contain the usual couplings such as the lens-lens correction and the correction to the Born 
approximation [101 1131 Hll HI US] The third line is the correction involved in the relation between the shear 
and the reduced shear, containing the coupling between the first-order convergence and shear [TO] . 

If we relax the small-angle approximation and we consider the full sky, there are other terms which become 
important T. The total geometrical contribution is 

^ _ ^(stan) I (corr) , (z) ^-,o^ 

/geom — /geom ^ /geom ^ /geom i V-^"/ 

where the standard terms 7goom'' are given in eq. (17), and the new terms, that we can consider as corrections 



^For comparison with a recent article in the literature on the subject, by rewriting in the first two lines of eq. \n\ 
the spin raising and lowering operators in terms of spatial gradients one recovers the second line of eq. (8) of [IT]. 
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to the small-angle approximation, are given by 







XsX 



2 rdx! + ^ ) vl/(x,f)^2^(x',x') (19) 



^0 



XsX XsX' 



2/ dx / dx'^V^*(X:2')^'*(x',2?') 
Jo Jo XsX 



Note that the second hne of eq. (17) contains only two transverse derivatives and its contribution to the 



shear is therefore of the same order as the corrections in eq. (19 1. Thus, one could be tempted to include this 
term in the geometrical corrections rather than in the standard contributions. The reason is that classifying 
terms by counting powers of $ and $ operators leads to ambiguities, as in the full sky those operators do not 



commute. It appears that the second line of ( 19 ) can be combined with the second term of the first line of 



eq. (17 1 to form a single term containing four spatial gradients, 

$$$^$^ + 2$m$m ^$0"^^$^ . (20) 

For this reason, we have decided to include this term in the standard contribution. 

Finally, the other contribution independent of the second-order components of the metric is the one due 
to the inhoniogeneity of the redshift of the sources. This induces a coupling between the photon redshift zs 
and the lens, 

^gcL = 1^-2 j^" dxi!{x,x) + ci,{xs,xs)-n-vs^ dx' 0^^{x\S') , (21) 

where H is the confornial Hubble parameter defined as "H = dlna/dr], vg is the peculiar velocity of the source 
and a dot denotes a derivative with respect to X: = ~ ^(^jj- Note that the inhomogeneity of the source 
redshift affects the convergence k already at first order in perturbation theory |39j . but its impact on the 
shear is only through a coupling with the lens, i.e. of second order. 

Since we are interested in ensemble averages over the random field </) and ipi we need to relate the velocity 
of the source vs to the primordial perturbation at first order. We choose to use the linearized Einstein 
equations to rewrite vs in terms of the Weyl potential at the source and its time derivative. In doing this 
we implicitly assume that the source is comoving with the dark matter. At large enough scales this is a 
legitimate assumption as all matter components are expected to move identically. However, at smaller scales 
this is clearly an approximation as galaxies typically belong to virialized halos and move independently from 
the dark matter field. At first order, the Oj-component of the Einstein equations reads 

^d,id^^ + n4>) = -%ixGp^v' . (22) 

Thus, using that ' = —9,,, the dark matter velocity at the source reads, at first order, 



^ LI'S V V— -r ^r,o^ 
TT2^ , (23) 



2 asV(-^s+Hs0s) 

where Q,m is the critical density of matter today, 0,^ = 87rGpm^o/(3^^o )> ^^^'^ ^® have used that p,„ = Pmfl/a' 



Note that the lens-lens correction and the correction to the Born approximation in the first line of eq. ( 17) 



can be written as the derivative of a deflection angle. The two shear components associated with these terms 



6 



are consequently related to the convergence and rotation part of the magnification matrix that contains only 
two degrees of freedom [4^ . On the other hand, the geometrical and redshift corrections cannot be written as 
the derivative of a deflection angle. Hence, due to these relativistic effects the magnification matrix contains 
in general four degrees of freedom: the shear £'-modes are generically different from the convergence and the 
shear _B-niodes are different from the rotation. 



(X, nx)^ I j—^ gix) T{k) e^^^-">^ , (24) 



2.2.2 Dynamical couplings 

Before moving to the contribution coming from second-order metric perturbations, i.e. the dynamical contri- 
bution, let us specify the first-order potentials and the initial conditions. Since we will perform our calculations 
in a ACDM universe, the anisotropic stress vanishes and the traceless part of the ij-components of the Einstein 
equations imply that and are the same at first order, i.e. c/)^^-' = ip'^^^ — ^E*'^^. Thus the Weyl potential 
can be decomposed in Fourier space as 

(2^ 

where T{k) is the matter transfer function and g{x) is the so-called growth-suppression factor, defined as 
g{x) = D{x) / a{x)^ where D[x) is the linear growth function. Furthermore, $ is the primordial potential 
which represents the initial curvature perturbation generated during inflation. During matter dominance, on 
super-horizon scales it is simply proportional to the curvature perturbation on uniform density hypersurfaces, 
$ = — (3/5)C. For the minimal model of inflation C, is approximately Gaussian [IT] and in the following we 
are going to assume that <i> obeys perfectly Gaussian statistics. 

At second order we split both (j) and i]j into a Newtonian part equal for both potentials, denoted by 

(2) ... . (2) (2) 

4>N , and a relativistic part, denoted respectively by cf)j^ and "0^ • The Newtonian part dominates on small 
scales; on these scales the relativistic part is suppressed with respect to the Newtonian part by factors of order 
(Ha/k)^ and becomes relevant only on large scales. Thus, the Newtonian potential c/)^'' dominates in the small 
angle approximation and its contribution to the shear has been thoroughly studied in the literature [ini El 112] ■ 

(2) 

It is obtained from inserting the second-order Newtonian potential (j))^ into the first-order expression for the 
shear (16 1. This yields 

7Sr'= rdx^^^^0i^)(x,f). (25) 
Jo XsX 

('2') 

The Newtonian potential 0]^ has been derived in standard Eulerian perturbation theory [131 HI] and for 
ACDM it can be written as 

^^N\x,nx) = J ^0i^^(fc,x)e*"^ , (26) 

with 

'^^ = ' 3HiL ^^'i^^'"^ klklF.Mh, k..x)<^,,<^,, , (27) 

where the kernel F2,n is well approximated by 

1/1 , \ , ^1 ' ^2 ^1 , ^2 \ , 1 ^/^ r n2 



where e(x) — f [Pm(x)/ptot(x)]~^^^*^ [HI US] |^ Here and below we implicitly assume integration over repeated 

^Here we have neglected the contributions to F2,n from the radiation era studied in 46^. These can be easily added 
to the standard dynamical contribution without changing the other corrections. 
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momenta, i.e. J ^-^^3^(5_D(fc — A:i — fc2)- Note that the expression for F2,n above is exact in matter dominance, 

(2,) 

where e = 3/7 and F2.N is time independent. For ACDM the time dependence of (f>}^ cannot be factorized 
out in the growth-suppression factor g'^{x) in front of the integral but i^2,7V is mildly time dependent. The 
approximation above to characterize this time dependence reproduces very well the one given in 46 . 

(2) (2) 

On large scales the relativistic parts of the scalar potentials and become important. Also the 

vector and the tensor modes in the metric, respectively Wi and hij, generated at second order are relativistic 

(2) (2) 

corrections to the metric of the same order in aH/k as the relativistic potentials 0]j and 4'r ■ As they are 
intrinsically second order, they simply enter linearly in the expression for the shear, yielding 



- -^Hi^{x,x) + ihr{x,x)) - -2/l(Xs:^s) , 

where we have defined ujr = h^uJi and = n^Whij, the spin-1 part of Ui and hij respectively as lUi = e\u)i 
and ihj. = e^^n^hij, and the spin-2 part of hij as 2^ = e\e'_^hij. Here e+ = eg + ie^ where eg and are unit 
vectors orthogonal to the photon propagation h and to each other. The term containing 2/1 in the second line 
of this equation is a boundary term induced by the spin-2 part of the tensor modes, which account for the 



distortion of the shape at the source. Putting together the dynamical contributions of eqs. (25) and (29) we 
obtain [T] 

(stan) (corr) /on^ 

7dy„ = 7dy„ + 7dyn " (30) 



Let us give here the expressions for the relativistic components of the metric to be put in eq. ([29j). For 
a ACDM universe, these have been computed in j47j and their expressions are quite involved. However, as 
we will see in section [5] the effect of the dynamical relativistic terms on the shear bispectrum is extremely 
small. For this reason we can estimate their contribution by simply using the expressions for the relativistic 
components in matter dominance, where it is possible to factorize their time and momentum dependence. 
Indeed, in this case these metric components can be written, in Fourier space, as |2H 148) 



ct^fik.x) 



i^'i\k,x)^-^^'H'{k,x) , 



T{k^)T{k2) 

" /,(2)/ 



fci • ^2 — 3(fc • ki){k ■ k2) 



htj{k,x) 
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fc2 



[l-3ji(/cr7)/(M] 



r(fci)T(fc2) 

fc4 



- k' (fc2(fc • k2 

((fci • 



kl{k-ki) 



fel k-2 



ki ki){5'^ + k'y) 



^ k-y k'2 ^2 ^ k\ k'^ — 4 k'^ k^^ k\ ■ k'2 



(31) 
(32) 

(33) 
(34) 



In the expression for tensor modes, the spherical Bessel function ji (cc) is given by ji (a;) = sin(2;)/a;^ — cos(2;)/a;. 
One can check that this metric solution satisfies Gaussian primordial initial conditions. Indeed, on large 



scales and at non- linear order the gauge transformation from the C gauge [H] to Poisson gauge ( 14 ) is given by 
^ = —ip — |(/) in matter dominance. Thus, eq. (32) implies Gaussian initial conditions C^^-* = 0. Furthermore, 
we note that in the squeezed limit, i.e. when one of the two modes fci and fc2 is much larger than the other, 
the second-order relativistic components of the metric given above go to zero. This is another advantage of 



using the metric (14) with (j) and ■0 in the exponentials. 
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To summarize, the full-sky observed shear up to second order is given by the sum of five contributions, 
three geometrical and two dynamical [1]: 



^ = ^(>5tan) , (corr) , (z) , (stan) (corr) 
/ / gcom ' / gcom ' / gcom ' / dy n ' ' dy n 



(35) 



Their expressions are reported, respectively, in eqs. (17 1, (19 1, (21 1, (25) and (29). 



3 Angular power spectrum 

The aim of this section and of the following one is to give the explicit forms of the spectrum and bispectrum 
of the shear field. More specifically, we are interested in the spectra of the electric part of the shear field. This 
section focuses on the power spectrum calculation. To present the formalism it is convenient to introduce the 
lensing potential T through the relation, 

7(n) ^ $^T{h) , (36) 
where the shear correlation properties are entirely encoded into the statistical properties of T(n). 



3.1 From metric to shear spectrum 

To compute the lensing power spectrum we just need the shear at linear order in the Weyl potential. In this 



case, using the expression for the shear at first order, eq. (16), eq. (36) reduces to the usual definition of the 
lensing potential given in the literature (see e.g. 

rxs 



T(n) 



rfx^^^*(x,%) 

XsX 



(37) 



As we will see below when computing the bispectrum, this expression is inappropriate to describe the lensing 



potential at second order, while eq. (36) remains valid. 

The lensing potential T is a spin-0 operator. Thus, using the commutation relation (|6| one shows that 



(38) 



where to write the second line we have used that T is real at first order. From eqs. ([7| and ([s]), this equation 



implies that 20,im — -20,im and, using the definition of aE,im, eq. (12), that 



CLEM 



1(1,2) / dnY:^in)T{h) 



(39) 



Note that aB.im vanishes at first order (see eq. (13l). The angular power spectrum Cf" of the electric part of 
the shear, E(n), is defined from the 2-point function. 



(a 



E.lm'^EA'm' 



(40) 



As at this order the angular power spectrum of the magnetic part of the shear is zero, Cj^ — 0, there is no 
ambiguity in simply setting Ci = Cf in what follows. 



We then use the Fourier mode decomposition of eq. (24) for the Weyl potential and expand the plane 
wave in spherical harmonics, i.e. 



Jk-hx 



^'jiikx)YUn)Yr^{k) 



(41) 



Im 
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Combining eq. (41) with cq. (36) and (39), and integrating in dh with the orthogonahty relation of the 
spherical harmonics, we obtain an expression for aE.im as a function of the primordial potential <i>^, 



fliJ.im = K^: 2) / dxW{Xs,x) 
Jo 



where we have defined the window function W as 



x'x 



9{x) 



Finally, using eq. ( 42 ) in the definition of the C/ , eq. ( 40 ) , one obtains 

2k'^dk 



Ci = l\l,2) r^dxWixs.x) r^dx'Wixs.x') 
Jo Jo 



which can be rewritten in a more compact form as 

2k:^dk 



T'{k)P{k)ji{kx)ji{kx') , 



Ci = l\l,2) 



P{k) T\k) 



dxW{Xs,x)ii{kx) 



Here P{k) is the power spectrum for the primordial perturbation $, defined by 



{2nfd{k + k')P{k) 



For simplicity in our calculations we will only consider a scale invariant spectrum, 

P{k) = A^k'^ , 



(42) 



(43) 



(44) 



(45) 



(46) 



(47) 



although all the treatment can be easily extended to a non-zero tilt. In the following we will consider a 
ACDM cosmology with f2„i = 0.3 today and a Hubble parameter h — 0.65. Furthermore, to ease the 
numerical treatment, for T{k) we will use the BBKS fitting formula [50] . 



3.2 Numerical integrations: the Limber approximation 



In practice, integrating eq. (45) is numerically involved and time consuming. To compute the Ci it is common 



to employ the so-called Limber approximation [5T], which is valid for large I. This is based on the fact that 
the ordinary Bessel function J,y{x), related to the spherical Bessel function ji{x) by 



y = 1 + 1/2 



(48) 



grows monotonically from zero dX x = Q io x ^ v and then rapidly oscillates. For large I an integral of an 
arbitrary function multiplied by a Bessel function can be approximated by 



dx}{x)Mx) = .f{u) + 0{l/v'') , 
which can be written in the following form: 

^—f{k)ji{kx)ji{kx') = ^-^^^fi^lx) [1 + 0(1/^^^)] 

TT X 

Using this approximation one obtains the Limber-approximated Ci as 

r^s w\xs.x) 



Ci^i\l,2) dx 



X' 



T^{v/x)P{vlx) [l + 0{l/v^)\ 



(49) 



(50) 



(51) 
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Figure 1: Angular power spectrum computed using the exact expression eq. (45 I (solid black line), using the Limber 



approximation at first order as in eq. (51 1 {dotted- dashed blue line) and at second order as described in appendix [A} 
eq. ( |112[ ) (dashed red line). The redshift of the source is = 1. In the lower panel we have plotted the error made 
using the first- and second-order Limber approximation. The error made using the second-order approximation is 
smaller than 1% even at very small I. 



This expression approximates the exact C;, obtained by integrating the full expression in eq. (45), to better 
than 1% for / > 8. 

For small I this approximation fails. Since we are interested in studying the shear on the full sky and 



correlate fields lying at different x along the line of sight, we need to go beyond the Limber approximation (49 ). 
As shown in 52 , an integral of an arbitrary function multiplied by a Bessel function admits the series 
representation 

2 i/'' b 



dxf(x)J^(x) 



fix) 



+ 0(1/1.^). 



(52) 



This can be used to go at second order in the Limber approximation - i.e. to include corrections of order 
such as the second and third terms inside the bracket in the equation above - and to considerably improve 
the Limber approximation. The expression for the angular power spectrum at second order in the Limber 
approximation is given in appendix [Aj eq. (112). 

In fig. [l] we show the angular power spectrum computed using the exact expression eq. (45) compared with 
the Limber approximation at first and second order. By going at second order in the Limber approximation, 
one obtains an agreement with the exact expression to better than 1% for 1^2. 

In practice and for the calculations of the spectra and bispectra, except when otherwise stated, we will 
employ this improved approximation. 
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4 Angular bispectrum 



4.1 The £'-mode reduced bispectrum 



At second order, the lensing potential defined in eq. (36) is in general a complex quantity. Consequently 
2aim and -2a;m are not equal anymore. This implies that the second-order shear will generically have both 
i?-modes and i3-modes. 

As in the first-order case, we can express aE,im and as^im as a function of the lensing potential. From 
eqs. ([?]) and ^ and using the definition of aE,im and aB.im, respectively eqs. (12) and (13), one finds 



aE.ira^-l{l.-'2) j dnl^:„(n)Re[^2^(n)] = 2) j dn y4(n) Re[T(n)] , (53) 
aB,im^-l{l,-2) j dnYilXn)lm[$^-i{n)]^-l{l,2) j dhYi:^{h)lin[T{n)] . (54) 



These expressions relate the electric and magnetic part of the shear respectively to real and imaginary parts 
of the lensing potentials. 

The angular lensing bispectrum B^^^^^ is defined by 

where AT, Y and Z can be either the electric or magnetic parts of the shear, E and i?, and 

(56) 



h h h 
mi 1712 rn^ 



is the Wigner 3-j symbol. As B vanishes at first order, only expectation values containing three aE,im or two 
o-E.im and one aB.im will be non-zero. For parity reasons, the bispectrum involving three aEjm will be zero 
for h + h + h = odd, while the one involving two aE.im and one aB.im will be zero for + ^2 + ^3 = even. 

In the following we concentrate on computing the bispectrum of the electric part, which will be simply 
denoted as Bi^i^i^, defined as 

Bi,i,i, ■ (57) 

\mi 7712 7713 y 

Following the standard CMB notation in the literature [S3], we will concentrate on the reduced bispectrum 
bi^i^i^ defined as 

{aE,hmj_aE,l2m2aE,hm.3] = ^1^3 "hhh i \P°) 

where 



(59) 



Qhhh""' = / dhYl^r,^An)yhrn2{'^^)Yhn^An) 



{2h + l){2l2 + l){2h + l) (h h h\ I h h h 

mi m2 777,3 



477 \ 



is the Gaunt integral, which characterizes the angular dependence of the bispectrum. It naturally appears 
due to translational and rotational invariance on the sky, and enforces mi, m2, and m^ to sum up to zero, 
and li, I2 and ^3 to satisfy the triangle inequality. 
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To plot our results, it is convenient to introduce the weighted bispectrum, defined as 



ChCh + C'ljC/g + Ci^Ci^ 



(60) 



Note that the weighted bispectrum is independent of the power spectrum normalization A$ in eq. ( 47 1 



4.2 Geometrical contribution 

We start by computing the contribution of the geometrical terms to the bispectrum. In order to do this, let 



us rewrite eq. (53) by moving the operator 0^ in front of Yi*^{n) by integrating by parts twice, and by using 



eqs. (126) and (127) in appendix [Cj After these operations, the coefficient aE,im becomes 

aEMn = (-1)"+' / dh[.2K,ih) 7(n) + ^Yi^^in) 7*(n)] . 



(61) 



Now, let us initially consider only the standard geometrical contribution, eq. (17). Plugging this contri- 



bution in the above equation, replacing the Weyl potentials by their Fourier components using eq. ( 24 ) and 



employing the plane wave expansion (41 ), this can be rewritten as 



Iil2mim2 

( ^'Yilm, (n) $yU. (") + ^ ^K^n, in) $Yil^.^ (n)) + 2 $Yi\„,^ (h) $Y,\^^ (n) 
X dx rdx'W{Xs,X.)W{x,x')Mix,x') 



(62) 



+ 00Y*^^in)Y*^^{n) T'dx rfx'W^(x, , x)M^(x. , x')>^ (x, x') 
Jo Jo 



where is a function of x a-nd x' defined as 



M{x,x') = 



27r2 27r2 



T(/cl)r(fc2)jh(fclX)ij.(A2X')^/l™l(^l)^/2m,(fe)%%^ 



(63) 



In the fifth line [ ^ f^- ^ ] stands for the first four lines (on the right-hand side of the equality) after replacement 
of the spherical harmonic -2^;^ by 2^;^ and of all the spatial gradients by their complex conjugates. 

We have thus factorized out the projection operators from the time and momentum integrals. We can 
now make use of the properties of the spin- weighted spherical harmonics, described in appendix [Cj It is then 
straightforward to rewrite the integrals over the angle h in terms of the 3-j symbols using 



{2h + l){2l2 + l){2l3 + 1) h h h\ h h h 



An 



>Sl S2 S3/ \r71i 7712 



(64) 

Note that the 3-j symbol enforces si + S2 + S3 = 0. The fifth line of eq. (62), [0 ^ ], gives an identical 
expression as for the first four lines after sign change of the spin indices si, S2 and S3. Using the following 
useful property of the 3-j symbols, 



h h h 

.Si S2 S3 



(-1) 



h+h+h I '2 ^3 

\-Sl -32 -S3 



(65) 
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one finds that for h + I2 + I3 even this contribution is identical to the one of the first four hnes while for 
^1 + ^2 + '3 odd it is opposite so that all the contributions exactly cancel. Thus, as expected aE,im vanishes 
for li + I2 + I3 odd. 



The bispectrum is then obtained by correlating the second-order [aglom^]_E,im coeSicient in eq. (62) with 
the product of two first-order aE,im, whose expression is given in eq. (42). After taking the expectation value 
over the primordial perturbations using Wick's theorem and the definition of the primordial power spectrum, 



eq. (46), we obtain for the standard geometrical contribution, for li + I2 + I3 even, 

-1 



h h h 
, , 



3 -1 -2^ 

h h I3 
,0 2 -2, 



h h I3 



2 -2) 

+ i{hA)i{i2,mi + 2){h~i) 



h h I3 
, 1 1 -2, 



X3 



(66) 



X / dxidx2dx3 / dx'W{Xs,XiW{Xs^X2)W{xs,X3)W{x3,x')Cixi,X3;X2,x') 
Jo Jo 

-1 

+ l^{h,l)l{h,2)P{h,2) 



h h I3 
.0 2 -2, 



h h I3 
, , 



dxidX2dx3 / dx'W{xs,XiW{Xs,X2)W{Xs,X3)W{xs,x')C{xi,X3;X2,x') 



5 perms 



where 

^/ I II iii\ 

C[x,x ;x ,x ) 



2kldki 2kldk2 



T\k^)T\k2)ji, (fcix)jh {k,x')ji, {k2x")jh {k2x"')P{ki)P{k2) ■ (67) 



The bispectrum from the corrections to the geometrical contribution can be straightforwardly computed 



from eq. (19 1, analogously to the calculation for the standard contribution. For -I- Z2 -I- Z3 even we obtain 

fc^lu/.ia = i{h,mh.mh.2) dxidx2dx3W{xs,Xi)W{Xs,X2) 

Jo 

X W{Xs,X3)9{X3)C{Xl,X3;X2,X3) 

- 5(X3) + (Xs - X3).9'(X3)^ J dx'9{x')C{xi,X3;X2,x') 
I h h h \ I h h I3 \ 



Xs X3 

2l{h,2)f{l2,2) 



2 -2/ \0 



dxidx2dx3W{Xs^Xi)W{Xs.X2) ^^^^ 



L"'0 



\ X3 



+ ^'^^1 5(x')C(xi, X3; X2, x') 

XsX' / 



5(X3) 



Xs Jo 



rfx'w^(xs,x')C(xi,X3;x2,x') 



5 perms . 



The calculation of the bispectrum from the redshift correction follows from eq. (21 ) and proceeds as above. 



The only complication resides in the spatial gradient in the velocity term, eq. (23), projected along the line 
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of sight. Using h ■ Ve*'^ "^ = d^e^'''"^, this simply introduces a derivative with respect to x of the spherical 
Bessel function in the expansion of plane waves into spherical harmonics. In this case the bispectrum reads 

/ \ / \ ~^ 



h h h \ I h h h 



2 ~2J \0 OJ 

X / dxidx2dx3W{xs,Xi)W{xs,X2)g{x3) 
1 



nsxl 



dx'9'{x')C{xi,X3;X2,x') - 5(Xs)C(xi,X3;X2,Xs) 



(69) 



2as 



gixs) - 



9'{Xs) 
-Hs 



'D{Xl,X3\X2,Xs) 



where 



2?(x,x';x",x"') 



5 perms 



_ f 2kldki 2kldk2 



T\k^)THk2)]i, {k,x)ji, {kix')ji. (fc2x") 



djiAhx" 



-P{k,)P{k2). (70) 



ainx'" 

4.3 Dynamical contribution 

We compute here the bispectrum from the dynamical terms. Let us start by the Newtonian contribution, 



eq. (25 1. For this contribution the lensing potential is real and reads, using eq. (36), 



T 



(stall) 
dyn 



dx— -(t^N^ix^x) 

XsX 



(71) 



Plugging this expression with (27) into eq. (53), expanding the plane wave in spherical harmonics and inte- 
grating over the angle n using the orthogonality of the spherical harmonics we find 



3i?gf2„i JQ 



d%d% 
(27r)3 



c'xl^(Xs)X)5(x)a(x) 



d^ks 

-^YLih)ji{k3X) 



(72) 



<5(fc3 - fci - k2)T{k^)T{k2)^F2M{klM,X)^i,^u, 



where F2^n is given in eq. (28). 



Using eq. (42), the expectation value of three asjm reads 

{aE,hm,aEMrn,aE,l,m,)^iyT^ = - i {h , 2)i {h, 2)l{h, 2) 



3/720 



X / dxidx2dx3W{xs,Xi)W{Xs,X2)W{xs,X3)9{X3)a{x3) 

~2^ ^^-^'^ (fclXl)j(2 (^2X2)^(3 ik3X3)Yil^^ (fcl)^(2m2 (^2)1^;3™3 (fcs) 

xT\ki)T^ik2)'^F2,N{ki,k2,X3)Piki)Pik2){27:fSDih + k2 + h) 
5 perms . 

(73) 
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Now we can write the Dirac delta function in the last line of this expression as 

1 



(27r)3 



(74) 



X Yj',m;(^i)>1>;,(fc2)yi>!,(fc3) / dxx ji[{kix)ji'^ik2x)ji'Jhx) 



where we have used y = nx- Integrating eq. (73) over the angular directions d ki using the above equation 
and the orthogonality relations of the spherical harmonics, we find for the bispectrum 



dxidX2dx3W{xs , XiWixs , X2)W{xs > X3)ff (X3)a(X3) 



JL. 01-2 01-2 



2kfdki 2kidk2 2k^dk3 kfki^ 



TT TT 

,2„ 



T'{ki)T\k2)F2^N{ki, k2,k3;x3)Piki)P{k2) 



(75) 



dx X jh {kiXi)jh {kix)ji2 (^2X2) j/2 {k2x)jh (^3X3)^(3 (kax) 
+ 5 perms , 

where the kernel -P^.w is a function of the amplitudes of the three momenta ki, ^2, ^3, as easily shown using 



ki ■ k2 



{kl 



kf - kl)/2 in eq. (28). 



We now move to the dynamical corrections, eq. (29). For the rest of this section we will assume matter 
dominance, so that the growth suppression factor is g = 1 and the window function reduces to ^^(XjX') = 
(x ~ x')/(xx')- The first line of eq. (29) yields a real lensing potential. For the scalar relativistic contribution 
coming from the first two terms in eq. (29), i.e. {(jrj^' +i/)]^')/2, the computation is analogous to the standard 
scalar part above, the only difference being the kernel. The details of the calculation and the result can be 
found in appendix [B} 

The computation of the vector and tensor contributions is more involved since their kernels depend not 
only on the amplitude of ki and k2 but also on their directions. We present here the computation of the 
vector modes. The computation of the tensor modes, which is similar, can be found in appendix [B] The 



vector modes contain two types of terms. The first one comes from the first line of eq. (29). Using eq. (33) 
and after some manipulations it can be written, in Fourier space, as 



\^r{k:i) 



i2ai/2 T{ki)T{k2 



fel fe2 



ksikl + k2) — (kf k'^)'^ T' _ 2 



2fc| 



{k^ ■ h) — k2{ki ■ n) — k^{k2 ■ n) 



(76) 



The computation of the term proportional to k^ ■ n is analogous to the computation of the Doppler term in 
the redshift correction. Using that 

ik^i-fK^^^-^-^ ^d^e"^^-^^ , (77) 
the time derivative of Fourier modes translates into a time derivative of spherical Bessel functions in the 



expansion of eq. (74). The rest of the calculation follows the one for the scalar modes. 

The same trick cannot be used for the second and third terms proportional to ki-h, i — 1, 2, since it would 
require to split the exponential as e^'^^"^ = e^'^i'^e*'^^'^. This would introduce additional infinite summations 
over V and m' in the expansion of the exponential into spherical harmonics. The alternative solution is to 
expand directly the product ki ■ fi into spherical harmonics using 



(78) 



M 
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With this expansion we find that the aE,im coefficient from the second term in the brackets of eq. (76) reads 



X 



3 M=-l,o,l 



X YiMiki)Yi'^ni[{-ki)Yi'^m'^{-k2) / dxx Jz; (fcix)j;i (^2x)j/^ (fc3X)ji!, (^aXs) ■ 

^0 

When computing the 3-point function, one ends up with three spherical harmonics to integrate over fci , which 
give rise to another Gaunt integral. The summation over three Gaunt integrals can then be simplified by 
using that 

V / — , — 7n2 , — mjj ^ — ma ,mjj , Af mi , — A/ 

V 47r 47r \rni m2 m J 

^ \0 oj \0 oj \0 Oj\l's 1 l[j ' 

where < > denotes the 6-j Wigner matrix. The properties of this matrix impose that the only terms 

in the summation over l[ and that give non-zero contributions are those with = /i ± 1 and ^3 = Z3 ± 1. 
Using the same procedure for the third term in eq. ( 76 ) , we find that the bispectrum of the vector contribution 
— ■^uJr reads 

[b%7\hh ^^l{h,2)l{h,2)l{l,,2) dxidx2dx3W{xs.Xi)W{Xs.X2)W{xs.XzW\xz) 



dyn J'1'213 3ijQ 

2kldki 2kldk2 2fc|dfc3 2 



Vr TT TT 

{k\ + fcj)^! ~ (^1 ~ ^2)^ 



T\k^)T\k2)P{k^)P{k2) 



X / rfxx jiAkiXi)jiAkix).n2ik2X2)ji2{k2X)dx3ji3(hX3)jhik3X) 







i3=;3±i 

^Ul h l3]fll I'l l\fl3 1'3 Afl'l h I'sXfk h hV' 

^ [is 1 ^'ij V 0/ \0 oJ\o oy 

X / dxx'^jiAkiXi)ji[{kix)jiAk2X2)jh{k2X)ji'3ik3X3)ji'3ik3X) > 



5 perms 
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The second vector term comes from the second line of eq. ( 29 ) and reads 



^iw(nx) 



3i7o J (27r)3 (27r)3 

'"fc|(/cj + fc|) — (k\ — fc|)" 



2fc| 



(fca • e+) - fci(fc2 • e+) - kl{ki ■ e+) 



(82) 



Here the new terms are those proportional to k^ ■ 6+. To compute these terms we rewrite them as fca • e+ = 
— ^(fcs • n) and then we use the expansion in eq. (77). The operator $ acts on the spherical harmonics 
giving a spin-1 spherical harmonic iYi^n{n), that results in a Gaunt integral with spin 1. We also have terms 
proportional to ki ■ e+, i = 1,2. Again, we rewrite these terms as fc^ • e+ = — 0{ki ■ fi) and then we use the 
expansion in eq. (78 1. With this, we find for the bispectrum of the vector term — ^ 



r,(corr)i ^ 

I dyn r^jj 



X3 



Jq 

TT TT TT 

(fcl + fcj)^! ~ (^1 ~ ^2)^ 

2fc| 

X / dxx^Jh (^iXi)jii {kix)ii2 {k2X2)ji2 {k2X)jh ihX3)ji3 (^ax) 



4Z^,(/!, + l) 



h I'l 1 



^3 1 



0/ \ 0, 



I'l I2 ^3 

, , 



+ j2V^{V^ + l){li+l'^ + 2) 



ll I'l 1 \ / ^3 ^3 



0, 



1-1, 



I'l I2 ?3 

, , 



(83) 



X j dxx jiAkixi)ji[{hx)ji2{k2X2)ii2{k2x)ji'^{hx'i)ii'.^{hx) I 

+ 5 perms. 

The computation of the bispectrum from the tensor contribution is analogous to the one above and is 
reported in appendix |B) 
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5 Results 



We now evaluate the weighted bispectrum from the geometrical and the dynamical terms in a ACDM universe. 
Since we are mainly interested in the impact of the full-sky corrections, which become relevant when one of 
the modes becomes very long (of the order of the Hubble radius) , we compute the bispectrum in the squeezed 
limit, li I2 — h- We label the low multipole corresponding to the long wavelength by k, and the two large 
multipoles, corresponding to the short wavelengths, by Ig and Is + Als respectively (see fig. [2]). 




I2 = Is h — h 

Figure 2: Possible choice of Zi, I2 and Z3 in the squeezed limit. Because of the triangle inequality enforced by the Gaunt 
integral, we let Als vary from to Z;. 
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Table 1: Color coding the different contributions to the bispectrum showed in fig. jsj 

In fig. [3] we plot the various contributions to the weighted bispectrum for two different values of the 
small multipole li and two different values of the redshift of the sources Z5. The amplitudes of the different 
contributions depend on the redshift of the sources and on the value of the small multipole li and, as shown 
by this figure, there is a clear hierarchy between them. 

In the following we are going to discuss their amplitude and their shape, including their redshift depen- 
dence. Even though the bispectra shown in fig.[3]were computed using the second-order Limber approximation, 
for the following analysis we will make use of the first-order Limber approximation. Moreover, we will assume 
that the matter transfer function T{k) can be approximated by a power law, 

r(fc = (X f-l- V , (84) 



\X keq/ 

where k^q is the equality scale. Clearly, as shown in fig. ^ the power a depends on the scales that we consider 



At very small Vs, a — > for a scale invariant spectrum. Using this ansatz, we can employ eq. (51) to derive 
how the angular power spectrum Ci scales with I and with the distance of the observer to the source, Xs ■ We 
find 



I 



\XskeqJ 

The Limber approximated formulae for the bispectrum can be found in appendix [A] For a summary of the 
color coding in the figures and of the equations giving the expression of the various contributions see table [I] 
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Figure 3: The weighted bispectrum in the squeezed limit as a function of Is, for two different values of the small 
multipole and two different source-redshifts, i.e. zs = 1 (left-hand panels) and zg — 0.5 (right-hand panels). In the 
top panels we have plotted the weighted bispectrum for /; — 10 and < Als < 10 while in the bottom panels we have 
plotted the case for li — 4 and < Als < 4. The dynamical standard terms are in red, the geometrical standard terms 
in blue, the redshift contribution in magenta, the geometrical corrections in black and the dynamical corrections in 
green. For each term Als varies from (darker line) to its maximum (lighter line) by steps of 2. The sign in the legend 
corresponds to the sign in the large-Za limit. 



5.1 Geometrical correction 

We will start by describing the geometrical relativistic correction, eq. ( 68 1 , black line in fig. [s] This contains 
two terms that we study separately but will give identical behavior. The first term contains a total derivative 
^^(\l/(x, af)\E'(x', a;')) ^^'^ its Limber approximation is given by the first five lines of eq. (114). Using this 
equation we can find the scaling behavior of this term. 



^3 + 5 perms . (86) 

askeqj \XskeqJ 

Because the term above is a total gradient there is no modulation dependence on the angle between the short 
and long mode. We can compute the weighted bispectrum in the squeezed limit by using the expression for 



the Ci in eq. (85 1. Choosing li ~ li gives 



6,u.+A,, - ^^Ig^ - (1 + 0[l-^)) . (87) 

Thus, the weighted bispectrum for this term does not depend on the short mode Ig and gives rise to a plateau 
at large Z^, as shown in fig. [3j Furthermore, as it does not depend on the angle between the short and long 
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0.1 1 10 100 io3 10^ 

V(xfceq) 

Figure 4: The power a = dlnT/dlnk as a function of k/keq = l/{xkcq). For the cosmology chosen in this paper 
{h = 0.65 and Vim = 0.3) fee, = 42.7 Ho and x{z = 1) = 0.77 H^^ so that for I ~ 1000, a ~ -1.5. 



mode, it does not depend on A^^. However, its dependence on can be clearly seen by comparing the upper 
and lower panel in fig.jsj As expected from eq. (87), this contribution is large from smaller 



Note that the dependence on Xs drops out in the weighted bispectrum so that the latter depends very 
mildly on the redshift of the sources Z5, as can be seen by comparing the left- and right-hand panels in fig. [3] 
In order to better visualize the redshift dependence, in fig. [5] we have plotted the weighted bispectrum for a 
particular configuration [li — 10, Is — 1000 and < A/g < 10) as a function of 25. 

Let us discuss now the second term, containing \E'(x,a;) ^^^'(x', a;')j i-^- ^ coupling of a spin-0 with a 



spin-2 field. Its Limber approximation is given by the last five lines of eq. (114). By applying the same scaling 
argument as above we obtain 

T \ 2ai / , \ 2a2 



bhioh — -1 — h — 1 — Z2Cos(26l23) + 5 perms , 

\XskeqJ \XskeqJ 

where 6'23 is the angle between I2 and inside the triangle, and we have used the general expression 

1 

« (-l)^cos(s-023) , (89) 
valid for large /'s. If in the squeezed limit we take li = li and I2 — I3 = Is, the dominant term in the 




permutation is the first on the right-hand side of eq. (88 1 or the equivalent term with 2 and 3 exchanged. 
These are proportional to cos(26'23) - where ^23 is the angle between the two short wavemodes - which in the 
squeezed limit goes to 1. Thus, also for this term we expect no modulation on the angle between the short 



and long mode and for the weighted bispectrum we find the same scaling as for the first term, eq. (87 1. Again, 



the dependence on Xs drops out in the weighted bispectrum, whence it does not depend on the redshift. 
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5.2 Redshift correction 

Let us discuss now the redshift correction, eq. (69), in magenta in fig. [3] This contribution is of the form 
^{x,x) ^^^'(x', a?')(Q ^^'^ same arguments as for the second term above can be apphed (the first-order 
Limber approximation for the redshift correction is given in eq. (115)). Hence, the redshift correction is also 
constant in Is and is independent of Als, as shown in fig. [s] (magenta hues). However, the redshift dependence 
of this contribution is more complex than the one of the geometrical correction and we need to study it in 
more details. 



Indeed, as shown in eq. (21), the redshift correction is made of three contributions: a Sachs- Wolfe, an 



integrated Sachs- Wolfe and a Doppler term. Using eq. ( 69 ) we can infer the scaling with the source distance 
from the observer. In the squeezed limit the first two contributions scale as 

hij.+AL-ir' — ii + oii;')) , (90) 

where rj^, is the conformal time at the source. The extra factor rj^/xs with respect to the geometrical correction 
comes from the term in front of the first square bracket of eq. ( 69 ) . For the Doppler term one finds 

2 



7,U.+A/. (-) {l + 0{l-')). (91) 



Thus, for sources at very high redshifts, i.e. such that ry^ < Xs i the redshift correction is typically dominated 
by the Sachs- Wolfe and integrated Sachs- Wolfe term, and is subdominant with respect to the geometrical 
correction. For sources at low redshifts 77^ > Xs the Doppler term dominates over the Sachs- Wolfe 
and integrated Sachs- Wolfe terms, and the redshift correction is larger than the geometrical correction. The 
redshift of transition is around zg ^ 1.2 but its exact value depends on the particular configuration. This 
is shown by comparing the left- and right-hand panels in fig. |3] and more directly in fig. [5] where for each 
contribution we plot a particular configuration as a function of the redshift zs- 



*The Doppler term contains an additional gradient of '5, but this has no impact on the Wigner symbols and 
consequently on the Als dependence of this contribution. 
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A remark is in order here. Naively one would think that the contribution coming from the Doppler term 
vanishes. Indeed, in the small angle approximation the correlation between the Doppler term and the first- 
order shear is taken at the same time x ^ Xs and vanishes for two reasons: the first-order shear vanishes at 
the source because of the window function; the Doppler term is longitudinal to the line of sight while the shear 
is transverse, so that being orthogonal their correlation at the same time vanishes. However, the correlation 
of the lens with the source does not vanish immediately for x < Xsj but for a given mode I it remains large 
as long as Xs ~ X ^ x/^ after which it decays exponentially (see appendix of [2l])- In other words, the 
correlation between the source and the lens decays when the distance from the source is of the order of the 
typical wavelength. This also implies that the redshift correction is dominated by the configuration in which 



the Doppler term and one of the first-order shear terms in eq. ( 69 1 is in the long mode. 



5.3 Standard geometrical contribution 

We now study the standard geometrical contribution, given in eq. (66) and shown in blue in fig. [s] This 
contribution is given by the collection of the lens-lens coupling, the Born correction, and the non-linear 
correction induced by going from the shear to the reduced shear. All these terms contain four transverse 
derivatives of the Weyl potential and we expect that they dominate over the two contributions discussed 
above. 

Let us start discussing the correction induced by the reduced shear. This term is proportional to 
$'^'^{x,x) $ $'^{x',x') (see e.g. last line of eq. (17)) and hence describes the coupling between a spin-2 



and a spin-0 field. Using the first-order Limber approximation given in the last two lines of eq. (113), the 
bispectrum scales as 

'1 \ , I h 



biM, - ( h ] h cos(2023) + 5 perms , (92) 



where we have employed eq. (891 to rewrite the Wigncr symbols in terms of the cosine. If we chose li = Z 



the four dominant configurations in the squeezed limit are those containing /i, i.e., summing up all of them, 

I \ 2Q1 / I \ 2a2 

'1 \ , I h 



^hhh h —t- /2[cos(20i3) +cos(20i2) + 2cos(2023)] , (93) 

XXs'^eqJ XXs'^eqJ 

where we have used I2 — h- In the squeezed limit 023 — > and 0i2 — > tt — 613, so that we can write the scaling 
of the weighted bispectrum as 

kus+AL (1 + cos(20i,)) (1 + 0{l^^)) . (94) 

Thus, this contribution is /s-independent for large Z^'s but its shape depends on the angle between the long 
and the short mode. For A/j ^ Is, cosOis — Al^/li and 1 -I- cos(20/s) ~ 2{Als/li)^ . This can be seen in detail 
in fig. [6] (left panel) where we have plotted this contribution separately from the other geometrical standard 
contributions. One can check that in the large-l^ limit the weighted bispectrum vanishes for Alg = 0, and 
becomes positive when A^s goes to li. 

We can now concentrate on the lens-lens coupling and the correction to Born approximation, which are 
described by the first two lines of eq. ( |17| . They contain three types of couplings: couplings between a spin-3 
and a spin-(— 1) field, (^■^'^{x,x) 0'^{x\x'), couplings between two spin-1 fields, 0^'^{x,x) $'^{x! ix'), and 
couplings between a spin-2 and a spin-0 field, $^'^{x, x) $ $'^{x\ x'). These types of couplings generate three 
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Figure 6: The weighted bispectrum from the reduced shear contribution (left panel) and from the lens-lens coupling 
and the correction to the Born approximation (right panel) at zs ~ 1, in the squeezed limit: li — 10 and < Als < 10. 
Als varies from (dark blue) to 10 (light blue) by steps of 2. 



different Wigner symbols, written in eq. (66). In order to study the scaling of this contribution we can use 



the recurrence identities for the Wigner symbols to write 

/ h h h \ /(/3-2)(/3 + ; ' 
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(95) 



so that in the large I's limit we can rewrite all the Wigner symbols in terms of cosines using eq. (89). Thus, 



inserting these decompositions into the first-order Limber approximation of the bispectrum given by the first 

B 

'1 



seven lines of eq. (113), the lens-lens coupling and the correction to the Born approximation scale as 

2a 

bhhh ~ ( ] h ( \ h 



h 



(cos(36'i3) -I- cos(6ii3)) 



L h 



2 cos(26'i3) -f cos(26'i3) -f C0S(26'23) 

'2 



-5 perms . 
(96) 



Taking the squeezed limit of this expression is subtle: at lowest order in 1/^^ (i.e. at zeroth order in the 
expansion of the cosines in terms of the small angle between the two short wavemodes) the expression above 
vanishes. Thus, we need to compute this expression at least at an order ©(Zj^) higher. Taking into account 
all the permutations one finds 



Kuis+i^U - (1 - 6 cos^ 6i,) (1 + 0{l-'')) . 



(97) 



This contribution is shown in details in fig. [6] (right panel). Using that 1 — 6cos^ 6*/^ ~ 1 — GAl'^/lf one can 
check that it is positive for A^^ < 4. 

Summing the reduced shear contribution to the lens-lens coupling and the Born correction, we find that 
the standard geometrical contribution gives rise to a plateau as a function of Is. The rcdshift dependence of the 
standard geometrical contribution is similar to the one of the geometrical correction. However, its amplitude 
is much larger due to the absence of the If'^ suppression present in the geometrical correction. Moreover, it 
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depends on Alg, as shown in fig. [s] (blue lines). As shown in the lower- right panel of fig. [sj for very small k 
and low Z5, the redshift correction becomes of the same order as the standard geometric contribution. 



5.4 Dynamical contributions 

The largest contribution to the cosmic shear comes from what we have called the standard dynamical con- 
tribution. This is given in eq. (75 1 and shown in red in fig. [sj This contribution is proportional to the total 



transverse Laplacian of a scalar, 0'^(f)^^\xi^) (see eq. (p5|)), but the scalar is not a Gaussian variable: it 



^(2) 



is related to the Gaussian primordial perturbation $ by eqs. (27 1 and (28). We can compute the scaling of 
this contribution as we did above for the others. We obtain 



"hhh 



Xs^eq 



Xs ^eq 



I9 



where, from eq. (28) 



F2,N{h,l2) = ^il + e) + ■ 
l^ we take the squeezed limit by choosing li 

2ai 



F2,N{h-,h) + 5 perms 
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-I- — cos 6*12 + - (1 - e) cos 
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^Xs^eq J \Xs^Gq 

where we have used I2 — h- As in the squeezed limit cos^ia 
and cos ^13 cancel and we finally obtain 



li, there are four contributions that dominate, 

h [2F2,n{IiM) + 2F2,n{IiM)) ' 

-cos 6*12, the terms proportional to cos 6*12 



2 [(1 + e) + {l-e) cos' Oi,] (^) ' (1 + 0(1-')) 

\Xs / 



(101) 



Thus, as can be checked in fig. [3j the weighted bispectrum of the standard dynamical contribution does not 
depend on Is but depends on the angle between the short and the long mode. This contribution dominates 
over the standard geometrical contribution because it is enhanced by the term (?7o/Xs)^. This term decreases 
at higher redshift, as shown by comparing the left and right panels of fig. |3]and more clearly by fig. [5] 

The dynamical correction from scalar, vector and tensor modes, see eqs. (1191, (81), (121), ( [83| , ( |123 1 
and ( 125 ), is shown in green in fig. [sj left-top panel. For large Is it is at least one order of magnitude smaller 
than the geometrical correction. Only for Ig < 100 it becomes larger than the geometrical correction and the 
redshift contribution. We will neglect it in the analysis of the following section. 
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6 Comparison with primordial non-Gaussianity 



6.1 Local bispectrum 

Let us start by studying the effect that a primordial signal would induce on the lensing shear bispectrum. 
The most common type of non-Gaussianity considered in the literature is the local type, which is typically 
generated by multi-field models [Ml |55l EH |57] . For local non-Gaussianity, the first-order primordial potential 
$ (which in matter dominance is related to the primordial curvature perturbation ^ by $ = — (3/5)C) contains 
a primordial quadratic correction in position space |58j . 

^x) = ^aix) - f^°l^U^) , (102) 
where <I>g is a primordial Gaussian potential whose power spectrum is given by eqs. (|46]) and (|47| Q Using 



the above equation together with eqs. (46) and (47), we can rewrite the 3-point correlation function of $ in 
Fourier space as 

(%*fc.*fc3> = (27r)35(fci + k2 + h)B'°''{h,k2,k3) , (103) 
where the local bispectrum is given by 

S'-(fci,fc2,fc3) = -2f^lA% (-I3 + ^ + J-3) . (104) 



We can use eq. (42) to compute the angular lensing bispectrum generated by local non-Gaussianity. It is 



straightforward to show that this is given by 

^tlh^ - i{h,2)l{l2,2)lik,2) r'dxidx2dx3W{xs,Xi)W{xs,X2)W{xs,X3) T'^Xx' 

Jo Jo 

2kldki 2kldk2 2k^dk3 . ^, n ■ x • /. n ■ ^, f, f, (105) 
Jh(hx)Ji2[k2X)Ji3{hx)]h{kiXi)Ji2(k2X2)]h{k3X3) ^ ' 

TT TT TT 

xT{k,)T(k2)T{ks)B'°%k,,k2,k3) . 

Let us use the Limber approximation and the scaling arguments of the previous section to find the behavior 
of the local bispectrum as a function of the long and short modes. Employing again the relation T{k) oc I", 
one finds 

In the squeezed limit we obtain, for the weighted bispectrum, 

\ As '"'eg / 

Thus, on very large scales, where ai — ^ 0, the local bispectrum behaves as the geometrical correction. 
6.2 Contamination 

The estimator of local primordial non-Gaussianity is given by |60) 



^Note that we use the notation of [59] with a minus sign for /nl- This is related to /nl used in CMB computation 



(which is defined in radiation domination) by /nl = 10/9/] 
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Figure 7: Overlap between the local non-Gaussianity and the standard dynamical and geometrical contributions and 
the geometrical and redshift corrections. 



where Ai^/^ig is a combinatorial factor equal to 1 if the three Vs are different, to 2 if two of them are equal 
and to 6 if all of them are equal. If the local bispectrum B^J^i^ is computed for = 1, the normalization 
factor A'loc that makes the estimator unbiased {{£) = /^l) is 



- E A r r r ' (^^^^ 

Thus, if this estimator is applied to an arbitrary signal Bf^^^i^^, this will contaminate the measurement by a 
value of given by 

yloc \ ^ hhh hhh (110) 

'1'2'.3 

Before estimating the contamination to the local signal from each of the contributions studied in the 



previous sections, let us use eq. (1081 to define a natural quantity that characterizes the "superposition" 
between two different contributions. We define the cosine between two bispectra as [ST] 

cos(r,X) ^ } J2 ' (111) 



where the factors Nx,y are defined analogously to A^ioc in eq. ( 109 1. This quantity ranges from —1 to 1 and its 
absolute value tells us how much two signals are orthogonal, cos(y, X) ~ 0, or superposed, cos(Y,X) ~ ±1. 
The cosine between the local non-Gaussianity and the contributions to the lensing shear from second-order 
perturbations is plotted in fig. [7] as a function of the redshift. For the calculation, we have summed over 
all Vs (with even sum) from 2 to Zmax = 1000. As expected, the geometrical relativistic corrections and 
the correction due to the redshift perturbation overlap very much with a local non-Gaussianity, contributing 
with a negative Indeed, in the squeezed limit both these corrections have the same I behavior and no 

dependence on the angle between the short and long mode. The overlap with the standard contributions is 
much milder. This is due to the fact that in the equilateral limit its amplitude decays less rapidly than the 
one of the local shape and because part of the standard contributions depend on the angle between the short 
and long mode, and this part averages to zero when convolved with the local signal. 
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Figure 8: The contamination to a primordial local non-Gaussian signal due to the standard dynamical and geometrical 
contributions and the geometrical and redshift corrections. 



We can now use eq. (110) to compute the non-Gaussian contamination to the local signal from each of 
these four contributions. The that an experiment looking for a local signal would measure from these 
contributions is shown in fig. |8]as a function of the redshift of the source zs- As expected, the geometrical 
correction leads to an f^f^ ^ 1 and roughly constant with zs ■ Even though the redshift correction has the same 
I dependence as the local signal, at low redshift it is enhanced by a factor {rjg/Xs)'^ ^-nd the contamination 
to f^l^ becomes surprisingly large, > 10 for zs < 0.3. The overlap between the standard geometrical 
correction and the local signal strongly depends on the redshift of the sources, and the of contamination 
inherits the same behavior. In this case /^l ^ 100 for zs > 0.8. Finally, as the standard dynamical 
contribution gets larger and larger at low redshift, its contamination to the local signal does the same. For 
zs < 1.7 we have f^l > 10". 
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7 Conclusion 



In this paper we present an exhaustive calculation of the full-sky bispectrum of the shear field as it would be 
built by an observer measuring shapes of galaxies at a given redshift. The calculation is performed at tree 
order (i.e. neglecting loop corrections), assuming that the initial curvature perturbation <^ obeys Gaussian 
statistics, such as after single-field inflation. Furthermore, we study the relative importance of the different 
couplings and we compare them with primordial non-Gaussianities of the local type. 

Our derivation fully exploits and completes the recent study of Ref. [1], where the full-sky second-order 
expression of the shear field has been computed. This article derives all second-order terms in the metric 
fluctuations that contribute to the shear field and it is thus the natural starting point of the computation of 
the tree order bispectrum presented here. In this paper we focus our calculations on the bispectrum of the 
scalar, i.e. electric, part of the shear field alone. This is somewhat restrictive since the second-order shear 
field contains a pseudo-scalar, i.e. magnetic, mode. We leave the study of the i?-mode and of the correlators 
where it is relevant for the future. 

The computation of the bispectrum from the formal second-order shear expression is an involved exercise, 
which includes two steps. The first consists in the formal derivation of the bispectrum in harmonic space. 
The electric part of the shear field itself is defined from its harmonic space expression, in a way similar to 
the mode decompositions of the CMB polarization. Its bispectrum is then computed with the help of the 
decomposition properties of the spherical harmonics and their spin weighted extensions. These results are 
presented in Sect. |4]of the paper where the contribution of each term is fully computed. Then, the second 
part consists in numerically evaluating the terms that we have found, in the context of the ACDM model. 

Two approximations are commonly employed for weak lensing calculations, namely the Limber approx- 
imation [52], which consists in neglecting longitudinal modes, as they tend to average out along the line of 
sight, and the flat sky approximation, which consists in replacing a small portion of the spherical sky with 
a 2-dimensional planej^ In order to be consistent with our full-sky treatment, it is a priori not possible to 
employ neither of these two approximations. Abandoning the Limber approximation makes the computations 
numerically challenging due to the large number of integrations involved and, as a consequence, only a limited 
number of configurations could be exactly computed. In order to circumvent this problem we have used the 
higher-order Limber approximation discussed in jSlj . Apart from for a specific type of contributions, we 
have found that sufficiently high accuracy can be reached with the help of the second-order approximation. 
The procedure that we have used to implement this approximation to the angular power spectrum and the 
bispectrum is presented in appendix [A] 

Our results are described in Sect. 5. In this section we give a precise account of the relative importance of 
the various contributing terms, as well as of their dependence with the multipole configuration. We distinguish 
the different contributing terms as geometrical and dynamical, depending on whether they originate from pure 
general relativistic effects on the line-of-sights or due to the dynamical evolution of the metric field. Note 
however that this distinction is somewhat arbitrary, since it depends on the gauge choice. Furthermore, in each 
case we distinguish the standard terms as those which survive at small angular scale, from the non-standard 
ones that are mainly relevant at large angular scales. 

We discuss the contribution to the bispectrum of the different nonlinear couplings focussing on the squeezed 
limit configuration, where the non-standard general relativistic corrections are more likely to play a significant 

^Note that these two approximations are not equivalent and one does not necessarily imply the other. 
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role. In Sect. [5] we give a precise account of the relative importance of the various contributing terms as well 
as of their dependence with the multipole configuration. We find that the general relativistic effects, i.e. the 
non-standard coupling terms, are typically smaller than the standard non-linear couplings. However their 
relative importance increases at smaller source redshifts. Among the new couplings, the dynamical couplings 
are always much smaller than the others, while those due to the inhomogeneity in the redshift of the source 
dominate. The latter can even become of the same order of the standard couplings when the redshift of the 
sources is below 0.5. These results are summarized in fig. |3]and[5| 

In Sect. |6] we finally compute the corresponding level of contamination induced by the standard and 
non-standard terms on the amplitude of primordial non-Gaussianities of the local type. Note that, as in the 
squeezed limit the standard non-linear couplings depend on the angle between the short and long mode, they 
differ significantly from a local signal. For 25 > 0.8, the contamination to f^{] of the standard geometrical 
couplings amounts to f^l > 100. On the other hand, the relativistic corrections induce a non-Gaussianity 
which is mainly of the local type and can contaminate the search for a primordial signal by f^^^ '-^ 1 to 
/nl ^ 10 when the source redshifts vary from 1 to 0.3. 

Our investigations did not go as far as computing the signal-to-noise ratio for the bispectrum induced by 
the general relativistic corrections. However, given the much larger amplitude of the standard contributions, 
we expect it to be rather small. This might not be the case when, besides the only iJ-mode three-point 
correlators we have investigated here, quantities involving i?-modes are also included. Relativistic corrections 
indeed generically induce B-modes with an amplitude a priori comparable to the iJ- modes [1]. This is at 
variance with the standard dynamical couplings that induce second-order £^-modes only. It ensures that the 
correlations between E and B modes are only non vanishing in the presence of general relativistic corrections 
making such quantities more likely to be observable. We leave such calculations for a future study. 
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A Limber approximation 

A.l Power spectrum in second-order Limber approximation 

The expression for the power spectrum with the second-order Limber approximation can be straightforwardly 



computed by applying the series expansion given by eq. (52 1 to the expression for the angular power spectrum, 



eq. ( 45 1 . The result is 



dXrr.2 
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T\v/x)P{i^lx)W{x) 



w{x) , x'W'ix) + -,x'w"'{x) 



(112) 



where W(x) = W{xstX)I \/X- By dropping the l/v^ term, one easily recovers the first-order Limber- 
approximated result given by eq. (51). 
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A. 2 Geometrical terms in first-order Limber approximation 



We present here the solutions for the bispectrum at first order in Limber approximation. In our numerical 
code we go beyond this approximation and we use the second-order Limber approximation. However, the first- 
order approximation is suflacient to understand qualitatively the scaling of the bispectrum in the squeezed 
limit. 

The bispectrum from the standard geometrical terms reads at first order in Limber approximation 
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where Ui = li + 1/2. 
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The bispectrum from the geometrical corrections is at first order in Limber approximation 
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The bispectrum from the redshift corrections is at first order in Limber approximation 
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Note that for the redshift corrections we cannot use the Limber approximation on the spherical Bessel functions 
ji2 ( vXs) and (^"Xs)- Indeed in this approximation x! — Xs and the bispectrum vanishes. Hence for 



these terms we perform the exact integral on x! ■ 
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A. 3 Standard dynamical term in first-order Limber approximation 



By applying the first-order Limber approximation given by eq. ( 50 1 to eq. ( 75 ) , it is straightforward to compute 



the Limber-approximated expression of the standard dynamical term. 
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B Dynamical contributions 



We present here the computation of the dynamical scalar and tensor terms. The vector terms are shown in 
section 14.31 



Let us start by the scalar terms. The sum of the two relativistic potentials gives, using eqs. (31) and (32) 

\ (0^''(fc3)+V4'^(^3)) =nfcimfc2)F2,s(^l,fc2,fc3)<ffe^$fe, , 
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For the bispectrum one finds 
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We now compute the tensor terms. We start with the contribution coming from the first line of eq. (|29|). 

ons one fi: 

T{k,)T{k2 



Using eq. (34), after few manipulations one finds 
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The computation of this contribution is similar to the one of the vector contribution, cq. (76 1. The only new 



term is the one proportional to the product (ki ■ n){k2 ■ fi). We use the expansion of eq. ( 78 ) twice. As a result 



in the computation of the aE,im we obtain an integral over d'^k^ with four spherical harmonics. This integral 
gives rise to a product of two Gaunt integrals. We then compute the bispectrum regrouping the appropriate 



Gaunt integrals such that the summation in eq. (80) can be applied twice. This generates two Wigner 6-j 
symbols, with terms of the form li ± 2. With this the bispectrum from the tensor contribution — |/ir-r is 
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+ 5 perms . 



We then compute the tensor term in the second hne of eq. ( 29 1 that reads 
10 f d% d^kid^k2 
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(121) 



(122) 



The computation of this term is similar to the one of the vector term in eq. ( 82 ) and we find for the 
bispectrum 
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Finally we compute the bispectrum of the boundary tensor term, i.e. the last term in eq. (291 
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+ 2(4 • e+)(A:^(fci • e+) + fc?(fc2 • e+)) - 2/c^(fci • e+)(4 • e+) 
The computation of this term is similar to the previous ones and its bispectrum reads 
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C Spin- weighted spherical harmonics 

The spin- weighted spherical harmonics sYim{n) are related to the standard spherical harmonics Yim{fi) by 

l{l,s),Yim{h)=0'Yim{n) , (126) 

s)^sYi„,{n) = i-iyfYi^{n) , (127) 
for s > and \s\ < I. The coefScient i{l,s) is defined as 



'('••) -Vl^^ (128) 

Let us recall here a few useful relations obeyed by the spin-weighted spherical harmonics, 

sYLin) = (-l)"+^_,y,,_„(n) , (129) 
$$Yi^{n) = l)Yi^ = -l{l + l)Yim{n) , (130) 

j dh sYim{n) sYi*^,{h) = 5w5mm' , (131) 

E sYim{n)sYil,in') = S{h - n') . (132) 

Im 

The operators and satisfy the commutation relation 

(0$ - $$) sYl^ = 2s ,Yira , (133) 

such that 

a^f V. - , 



fY,rn = -§ri^0Yirn. (134) 



The same relation holds when and are interchanged. From this relation we deduce that 

0''0^ Yim = l''{l,2)Yi^. (135) 

Dropping the boundary term when using integration by parts is valid as long as the integrand has spin zero, 
i.e. / dh 0{[spm = -1]) = J dn 0{[spm = 1]) = 0. 
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